Preferred sizes and ordering in surface nanobubble populations 
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Abstract 

Two types of homogeneous surface nanobubble populations, created by different means, are 
analyzed statistically on both their sizes and spatial positions. In the first type (created by droplet- 
deposition, case A) the bubble size R is found to be distributed according to a generalized gamma 
law with a preferred radius R* = 20 nm. The radial distribution function shows a preferred spacing 
at ~ 5.5R*. These characteristics do not show up in comparable Monte-Carlo simulations of 
random packings of hard disks with the same size distribution and the same density, suggesting a 
structuring effect in the nanobubble formation process. The nanobubble size distribution of the 
second population type (created by ethanol-water exchange, case B) is a mixture of two clearly 
separated distributions, hence, with two preferred radii. The local ordering is less significant, due 
to the looser packing of the nanobubbles. 
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The first atomic force microscopy (AFM) observations of spherical cap-like soft domains 



at the solid-liquid interface |, later termed "surface nanobubbles" , identified two 

typical, yet poorly understood, nanobubble characteristics: long-term stability and huge 
nanoscopic contact angles (on the water side). Later experiments confirmed these puzzling 
features of surface nanobubbles, and focused on verifying their gaseous nature by correlating 



the nanobubble densities with the gas concentration in the liquid 



8|. Recently, the 



gas content of the bubbles was identified explicitly by infrared spectroscopy measurements 
in combination with AFM 0, Q. Other studies investigated the effect of surface active 
H |, salts [ll|, substrate morphology 12], or electrolysis K^Q on the appearance, 



solutes 



stability, and shape of surface nanobubbles. While the number of experiments support ing 



the notion that the observed structures are indeed surface nanobubbles, is increasing 



16 



13, 



18 



19 



20 



21 



221 ]. no consensus has been reached concerning the mechanism which 



stabilizes the bubbles (see 



231 ] and references therein). One of the hypotheses, recently put 



forward in 23j, is based on a non-stationary equilibrium between a gas outflux (through 
the gas-liquid interface) and a gas influx (at the three-phase contact line), and predicts a 
preferred nanobubble radius as a function of gas concentration and contact angle. 



In this paper we want to test the prediction [23j of a preferred radius R* and its depen- 
dence on the gas concentration. Our good bubble statistics allow us to extract statistical 
properties of the whole nanobubble population. The analysis shows not only a preferred ra- 
dius, but also a preferred spacing between the bubbles, suggesting a structuring mechanism 
between individual bubbles. 

As substrates small pieces diced from a Si(100) wafer are used, which are subsequently 
cleaned, coated with a monolayer of lif ,liy,2/7,2/7-perfluorodecyldimethylchlorosilane and 
analyzed following the procedure described previously []jj]. The static macroscopic contact 
angle is typically around 92°. The substrates are then mounted in an atomic force mi- 
croscope (VEECO/Digital Instruments (DI) multimode) equipped with a NanoScope Ilia 
controller (DI, Santa Barbara, CA) and measured in tapping mode in water using a DI liquid 
cell and V-shaped SisN4 cantilevers, with spring constants of 0.3 — 0.5 N/m (Nanoprobes, 
DI). The amplitude setpoint was chosen as high as possible, typically > 90%. The size 
of the nanobubbles is extracted from the raw AFM topography images by application of 



a height-threshold 



141 ]. which yields the location and radius R of each nanobubble. The 



results are corrected for the finite size of the tip (R t i P =20nm), as done elsewhere [21]. We 
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FIG. 1: AFM topography images of the solid-liquid interface of the substrates. In a) and b) 
gas-equilibrated MilliQ-water was put on the substrate without explicit use of local oversaturation 
(case Al and A2, resp.). In c) and d) the result is shown after a local and temporal oversaturation 
has been applied (case Bl and B2, resp.). Each scale bar corresponds to lfim. 



note that the tip correction does not affect the conclusions of this paper qualitatively. 

The populations of surface nanobubbles are created in two different ways: in case A, 
a drop of gas-equilibrated Milli-Q water is put on the substrate, while in case B a finite, 
temporal local gas oversaturation (by flushing ethanol away with water [l|, [ill, [3]) is em- 
ployed to explicitly stimulate nanobubble formation. In both cases two typical images are 
selected which were suitable for further statistical analysis (see Fig. [1]). In Case A (Fig. [T^i 
and b) a dense coverage of relatively small and rather uniformly sized nanobubbles is ob- 
served. This observation is not evident as not all labs find the 'spontaneous' occurrence of 
nanobubbles (see for instance the remark in [lfj and references therein). Only incidently, 
some larger nanobubbles are visible, which are present next to a bubble- free area. Presum- 
ably, smaller nanobubbles have merged to these larger objects. A mixed population of both 
small and large nanobubbles can be created when a forced local oversaturation is applied 
temporally [uj, as shown in Fig. [Ifc and d (cases Bl and B2, resp.). After the local gas- 
oversaturation, the bulk gas concentration is restored to its equilibrium value. In addition, 
the bubbles have been exposed to a single shockwave, as described in [19(] . We noticed that 
the large nanobubbles did not vanish or shrink during the course of the experiment (i.e. 
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within a few hours). 

The experimental probability size distributions P(D) present in case A and B are shown 
in Fig. [2J The bubble sizes clearly show a maximum at a particular diameter value, which 
we denote as the preferred diameter D*(= 2R*). In case B there are even two peaks, 
corresponding to two preferred radii. To obtain the value of D* the experimental size 
histograms were fitted with a generalized gamma distribution (GG) in case A and with 
a mixture of a GG and a Gaussian distribution in case B. In Case A the GG distribution 
which best fits the experimental results (Fig. [2]) is 

Pa(D) = ^0^ De M-(PAD) 9 *h (1) 

where is Euler's gamma function, and (3a and 9a are shape parameters which are fitted, 
yielding (3 A = (1.73 ± 0.07) • lO^nm -1 and 6 A = 2.37±0.17. As the value of the exponent of 
D in front of the exponential was found to be very close to 1, it was fixed to 1. The maximum 
(or the mode) of Pa(D) is formed at D* A = 1/ ' (5a0 1 J l ° a = 40 ± 2nm and the mean diameter 
(D)a = 47 ± 2 nm. The standard deviation of the size distribution is a a = 23 ± 2 nm. 

In case B the total probability distribution could be fitted with a mixture of a GG 
distribution with the same form as that of Pa{D) and of a Gaussian distribution: 



9 bP% n _ r , a m * H1 l-« / (D-D 



\2 



Pb ^ = "F7^^ ex P{-(/W s } + ^exp - V . (2) 

As in case A, the exponent of D was found to be close to 1 and fixed to that value. The 
fitted parameters are then a = 0.69 ± 0.04, f3 B = (1.20 ± 0.07) • lO^nm -1 , B = 2.8 ± 0.3, 
D%2 = (D)b,2 — 224 ± 9nm, and &b,2 = 48 ± 7nm. The characteristics of the small 
nanobubbles are then D* B 1 = 58 ± 4nm, (D)b,i = 64 ± 4nm and <jb,i = 29 ± 4nm. 

The observation of two co-existing but clearly separated sets of bubbles has not been 
reported or predicted before. The larger nanobubbles are created during the temporal 
gas oversaturation in the water during the exchange process, in agreement with previous 



observations H], while we hypothesize that the smaller ones are formed once the saturated 
conditions are restored. Notice that the smaller set of bubbles in case B is fairly similar to 
the population in case A in both the shape of P{D), as well as the order of magnitude of the 
respective maxima (40 ± 2nm and 58 ± 4nm,resp.). Remarkably, these maxima are close 
to the experimental result of Simonsen et al. [17], who found a normal distribution of sizes 



4 



20 40 60 80 100 120 140 




20 40 60 80 100 120 140 

D (nm) 

50 100 150 200 250 300 350 400 





D (nm) 

FIG. 2: Probability distribution of the nanobubble diameter D in both case A (top) and case B 
(bottom). Each case is represented by two unique images (1 and 2, resp.), of which the total size 
distribution is shown. The bars depict experimental data, the lines show the best-fitted probability 
distribution. 
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FIG. 3: Positions of the nanobubbles in case A2 (left) vs. Monte-Carlo simulations of a random 
packing of hard-disks with the same size distribution and density as experiment A2 (right). The 
experimental positions show much more structure than the simulated bubbles. 

with D* = 66 nm under identical lab conditions (i.e. gas-equilibrated Milli-Q water put on 
surfaces with a static contact angle 90°). 

The homogeneity of the nanobubble coverage depicted in Fig. CD suggests local structuring 



5 



0.5 




2 4 6 8 10 

r/<R> 



0.5 




2 4 6 8 10 

r/<R> 



0.5 



0.5 




2 4 6 8 10 

r/<R> 




2 4 6 8 10 

r/<R> 



FIG. 4: Radial distribution functions g(r) as a function of r normalized by the mean radius (R) 
for case A (top) and case B (bottom). Blue line: experiment; red line: Monte-Carlo simulations 
of a random packing of hard-disks with the same size distribution and density as the associated 
experiments; dashed line: Poisson point process; dash-dotted line: determinantal point process. 

of the bubbles. To test this idea quantitatively, Monte-Carlo (MC) simulations of a random 
packing of hard-disks with the same size distribution and density as in the experiments 
are employed. For case A2, the nanobubble center positions in both experiment and MC- 
simulation are depicted in Fig.[3j which shows that the experimental positions are much more 
structured than the simulated bubble positions. This effect is further shown by the radial 
distribution function g(r), which quantifies the probability of finding a bubble at a radial 
distance r from another bubble, and the nearest neighbor distribution function D]yN(r), 
which gives the probability of finding a nearest neighbor of a nanobubble at a distance 
less than or equal to r [25|. The plots of g{r) and D NN {r) are depicted in Fig. [4] and 
Fig. El respectively, for both the experimental and MC-simulated positions. In addition, the 
figures show the distributions for a Poisson point process (where neither steric nor repulsive 
interaction is p resent), and a determinantal point process with a very soft repulsion between 



the points 



For case Al and A2 the experimental curves in Fig. H] show a significant peak in g(r) 
at r ~ 5(R) ~ 5.5-R* while, interestingly, this peak is absent in the corresponding MC- 
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FIG. 5: Nearest neighbor distributions Dj\u^(r) for the four cases A1-B2. The legend is the same 
as in Fig. [H 

simulations and the determinantal point process. This shows that there is a preferred spacing 
between the bubbles present in both case Al and A2, which is not only steric and stronger 
than the 'soft' repulsion represented by the determinantal point process. The regularity of 
the bubble positions in case Al and A2 is also shown in the plots of D^^lr) (Fig. the 
experimental curves are on the right-hand side of the MC-curves. Notice that the D^^ir) 
curves of the MC simulated cases Al and A2 (red lines) are close to those given by the 
determinantal point process (dash-dotted lines) although in the MC simulations hard disks 
are used without any mutual interaction. The similarity is not seen when the MC simulation 
utilizes a single disk size. Hence, the effect of the size distribution looks like an effective soft 
repulsion. 

In contrast, for cases Bl and B2 no significant difference is observed between experiments 
and MC-simulations, in both g(r) and D^^ir). The reason could be that the statistics is 
too poor: case A counted three times as many bubbles B. Another, more likely 

reason could be that the number densities in case B are too low for structuring effects to be 
present. In case Bl and B2 the number density was 13.8 per /mi 2 on average, while in case 
Al and A2 this was 70.7 per /zm 2 , more than a factor of five difference. 

In summary, it is demonstrated for two types of surface nanobubble populations that 
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nanobubbles (i) show a preference in size, and (ii) show a preference in spacing. The first 
observation shows up in both cases, while the second observation only shows up when the 
number densities are large enough. In case A the size distribution is found to be distributed 
according to a generalized gamma law. A very similar size distribution is present in case 
B, where in addition a larger set of normal-distributed nanobubbles is present, which were 
created most likely during the temporal gas-oversaturation in the water. These findings are 
consistent with the hypothesis of a uniform stabilizing mechanism leading to a preferred 



radius, as put forward in 



231 ] . Comparisons with MC simulations show that densely packed 



nanobubbles do not reside randomly, but choose a position were it is easiest for them to 
be: away from each others vicinity. The physical mechanism responsible for this effect 
could be the limited availability of gas in the vicinity of an already formed nanobubble, 
prohibiting the nucleation of other nanobubbles nearby. Alternatively, nanobubbles could 
be formed instantaneously from the breakup of a homogeneous gas film into individual 
bubbles, analogous to the break-up of thin liquid films into surface patterns [261]. Thirdly, 
the ordering effect could result from a short-range repulsive force, e.g. due to surface charges. 
Further theoretical and experimental studies are required to unravel the precise formation 
mechanism of nanobubbles and their mutual interplay at the nanoscale. 
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